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Abstract 



I have investigated what coordinate systems, other than the Robertson- 
Walker (RW) system, there might be to describe a homogeneous and 
isotropic universe that are physical close to the origin, i.e., covariant 
transforms from RW that use physical clocks and rulers. I show that 
there are an infinite number of such systems, and that all such sys- 
tems must have the same variable photon velocity close to their origin, 
which is consistent with the RW system also having the same variable 
photon velocity. Implausible though it may seem, if the criteria for 
physicality are valid, it seems inescapable that this must represent the 
physical photon velocity c(t). The assumption of homogeneity makes 
it applicable to physical processes throughout the universe. I find that 
the analytic expression for it is c(t) = aH, the cosmic scale factor times 
the Hubble ratio in appropriate units. The existence of variability de- 
rives from the physicality criteria being applied to the transformed 
coordinates expanded from the origin out to the lowest order of the 
RW radial coordinate, and not to gravitational energy density, which 
determines only its magnitude. With the use of a generalized time for 
which the differential is c(t)dt, we obtain a generalized Lorentz trans- 
form which is invariant in a generalized Minkowski metric. When we 
use the coefficients of this metric in Einstein's field equation, we lo- 
cally obtain the generalized Minkowski metric. With this I find that 
the Hubble ratio, defined with this generalized time, has the same de- 
pendence on the cosmic scale factor as it does for a constant photon 
velocity. Some other "constants" that depend on the photon velocity 
must also change. 
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1 Introduction 

The cosmological assumption that the universe is homogeneous and isotropic 
in space docs not automatically lead to a constant photon velocity. In fact, 
the Robertson- Walker (RW) space-time line element [1] [2], which can be de- 
rived from the assumption of this symmetry, allows for a variable photon 
velocity c{t) which depends on the cosmic time t. Although there is strong 
historic precedent for assuming a constant photon velocity, it should be no 
more a surprise to us that it is variable than that the universe size is vari- 
able in time, a great surprise to Einstein when he discovered his equations 
led to it, which he tried to fix with his famous cosmological constant. In 
this paper I propose to show that the same assumptions lead to a variable 
physical photon velocity. I do this by showing that the infinite number of 
spherically symmetric physical coordinate systems all require the same vari- 
able photon velocity c{t). By physical I mean that they are obtained by 
a covariant transform from the RW coordinates and that their coordinates 
close to the spatial origin are measured by physical clocks and rulers. If we 
use a valid description of physicality, one would have to conclude that they 
must describe the physical photon velocity. Homogeneity is assured by re- 
quiring each coordinate system to be the same no matter where its origin is 
located in the universe. This means that the same variable photon velocity 
applies to all physical phenomena throughout the universe. 

We assume that RW system provides a good physical description of the 
universe and that its differentials can be measured by clocks and rulers. A 
covariantly transformed coordinates system also describes the same universe 
because by definition of covariance it satisfies the same field equation as the 
RW system. To determine whether the transformed coordinates are using 
physical clocks and rulers, we will use the Bernal criteria[8] to be assured 
that the differentials are as physical as the differentials of the RW system. 
We will require this condition to apply out to the lowest order of the radial 
RW distance when expanded away from the origin. Homogeneity requires 
invariance of the photon velocity under reflection: an outgoing photon from 
the origin will be an incoming photon to another close-in galactic point, 
which by homogeneity is required to equal the incoming velocity to the 
origin if measured at the same instant of cosmic time. We will see that this 
requires the metric of the physical coordinates to be diagonal close to the 
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origin. 

A variable photon velocity has far reaching imphcations. A large num- 
ber of other physical "constants" depend on it, so many of them presumably 
also change with cosmic time. In addition a variable c{t) no longer keeps the 
Lorentz transform invariant, so that the equivalency principle must be mod- 
ified. Following the suggestion of Magueijo [3], we introduce a generalized 
time t for which dt = c{t)dt in order that locally we have an invariant gener- 
alized Lorentz transform and a generalized Minkowski metric for use in the 
Einstein field equation. Of course, this generalized Lorentz transform will 
reduce to the normal Lorentz transform for events whose time separation is 
small compared to cosmic time in order that all the laws of atomic, nuclear, 
and particle physics remain valid. I will show that the defined physical co- 
ordinates obey the generalized Minkowski metric and generalized Lorentz 
transform close to the origin. In fact we could have used this metric as the 
starting physicality requirement instead of the Bernal criteria. 

We proceed by using the RW coordinates with c = c{t) applied to a ra- 
dial world line to find the radial 4-velocity (Sect 2) of a point on the frame 
of the covariantly transformed coordinate system in terms of both the RW 
coordinates and the transformed coordinates themselves. The requirements 
of physicality close to the origin lead to a variable c{t) = aE, the normal- 
ized cosmic scale factor times the normalized Hubble ratio (Sect 3). An 
infinite number of such coordinate systems do exist with these requirements 
(Sect 4). When we use a generalized time t we obtain a generalized Lorentz 
transform, a generalized Minkowski metric, and a generalized Equivalency 
Principle. This generalized time replacing the normal time in the field equa- 
tion locally produces the generalized Minkowski metric (Sect 5). This allows 
us to calculate a{t) and thus c(t) (Sect 6). The variability of other physical 
"constants" is discussed in Sect 7. 

2 The transformation from RW coordinates 

For a homogeneous and isotropic universe, by embedding a maximally sym- 
metric three-sphere (r, ^, ^ in spherical coordinates) in a four-dimensional 
space by including a fourth coordinate t, one can obtain a differential line 
element [4, page 412]: 

ds^ = citfdt^ - aitf[dx^ + r^du}\ (1) 
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where t can be taken as a universal cosmic time, and x is introduced as a 
co-moving curvilinear radial spatial variable (x > 0), for which 



{sinx, k = l, 
X, k = 0, 
sinhx, k = —1 



(2) 



A; is a spatial curvature determinant to indicate a closed, flat, or open uni- 
verse, resp., a{t) is a scale factor multiplying the embedded three-sphere, 
and du}^ = dO^ + sin^Odcj)'^. This is the RW metric generalized to a variable 
photon velocity c{t). 

It will be convenient to introduce the time related quantity t [3], which 
we will call a generalized cosmic time, defined by 



The alternative line element dr^ = ds'^/cP could also be used to define a 
covariant derivative. This would give the same c(t), but is less convenient 
to use. 

The RW coordinates define a differential ruler distance a,dx, but the full 
coordinate x is not measured by rulers, but by numbered galaxy positions. 
The proper coordinates T = t, R = ax are an attempt to supply a ruler 
measured distance, but these arc not physical in that the resulting metric 
does not give an invariant reficctivc photon velocity when expanded from 
the origin to first order in x- I have tried to find a transformed system which 
provides coordinates that use physical clocks and rulers and have a reflective 
photon velocity. This is not possible for all values of x (see Sect 4), but is 
possible close to the origin, i.e., out to a small but flnite value of x- (There 
are other coordinate systems, such as the full RW system, which describe 
physical phenomena whose coordinates do not use clocks and rulers, but 
my objective was not to find all such systems, but only those that do use 
physical clocks and rulers, which I have thereby called physical). 

We will take the origin of the transformed radius R to be the same as 
the origin of the RW radius: i? = Oatx = 0. As a first step in making 
the transformed system be physical close to the origin, we will require that 
the transformed time T is t at x = 0, since the time on clocks attached to 
every galactic point is t, including the origin. With no loss in generality, we 




(3) 



so that the line element becomes 



ds^ = dt 



a2(dx^ + rW). 



(4) 
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will assume that the angular coordinates Og and ^5 of the new transformed 
coordinates are given by the angular coordinates of the RW system. We will 
be initially confining our attention to only radial trajectories, so that the 
angular coordinates are fixed, and T and R will be functions only of t and 
X- Thus 

T = T{t,x), R = R{t,x) 

Os = e, 4>s = </>. ^ > 

The transform of the difi^erentials is given by 

dT = Ttdt + Tydx, dR = Rtdt + Rydx, 

(6) 

dOs = de = o, d^s = d0 = 0, 

where the subscripts (other than s) indicate partial derivatives with respect 
to the subscript variable. We will henceforth drop the subscript s. 

The general form for a metric expressed in transformed coordinates with 
spherical symmetry is given by [4, page 336] 

ds'^ = A^citfdT'^ - B^dR^ - 2Cc{t)dRdT - a^r'^duP, (7) 

where A, B, and C are implicit functions of (T,R) and explicit functions of 
{t,x)- We will use as the RW coordinates t,x,G,<P and as the transformed 
coordinates T,R,9,(f), where dT = cdT, so that the transformed metric 
becomes 

ds'^ = A^df^ - B'^dR^ - 2CdRdf - a^r^duj'^, (8) 
and the transform of differentials from {t,x) to {T,R) becomes 



dT = Ttdt + cTy^dx, 
dR= l;Rtdi+R^dx, 



(9) 



Let us consider a point on the transformed system. When measured 
from the RW system, it will be moving at a velocity given by 

v = «w(|)^, (10) 

To get the time component of the contravariant velocity vector = dt/ds, 
we divide eq 4 by dP to obtain 
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To get the radial component, we use the chain rule apphed to eqs 10 and 
11: 

dx ^ dx£ ^ 
ds di ds ca 

Since the contravariant velocity vector of a point on the transformed 
system in the transformed coordinates is (1/A, 0,0, 0), this will transform 

from the vector = ^ in the RW coordinates (7, ^,0,0) the same as 
df, dR in eq 9: 

= lR,j + lR^^V. ' 

where V = V/c. 

Manipulating the second line of eq 13 gives 

The photon velocity is given by ds = 0. For the radial velocity Vp in the 
transformed coordinates we obtain 



Since we want the photon velocity to have the same value in both direc- 
tions, C = 0. This makes the metric diagonal. Thus, we will consider only 

transforms that become diagonal close to the origin. (In this paper we will 
use the adjective "diagonal" to describe any characteristic of the coordinates 
with a diagonal metric). 
If we invert eq 9, we get 

di=^{R^df-cT^dR), 
dx = U--cRtdf + TtdR), 

where 

D = TtR^ - RtT^ = TtR^il + i>^), (17) 

using eq 14. With this in the RW metric (eq 4 ), we find the values of the 
coefficients of eq 8: 



1 - V 
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(1 + ^^)2 



(18) 
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cT-y 



and 



To make C = 0, 



C = -- 



TtR 



V 



aTt 



0. 



(19) 



(20) 



(21) 



This can be zero either because each of the two quantities separately ap 
proach zero, or because 

aTt' 



(22) 



We assume eq 22 is the vahd relation for physicality because if they sepa- 
rately approach zero without being equal, their sum will be zero only at the 
origin, but not a short but finite distance away. In Sect. 3 we show that V is 
linear in x with no second order term, so that the relation in 22 is expected 
to hold to at least the second order in x- With eq 22 the expression for A 
and B become: 

(23) 



and Vp becomes 



B = ^ 



±C^. 

aTt 



(24) 



(25) 



The two different equations for V, eqs 14 for distance and 22 for time, will 
provide the critical constraint on the expansion of T and R to lowest order 
of X that will determine the variable photon velocity (Sect 3). 



3 Diagonal coordinates which are physical close to 
the origin 

We assume that the RW differentials represent physical coordinates: physical 
clocks measure dt and rulers measure adx, and that the spatial origin of 
the two systems is the same, so that R{t, 0) = 0, and V{t, 0) = 0. We 
use the criteria developed by Bernal et al[8] to determine whether the two 
coordinate systems use the same clocks and rulers. In their notation for xi 
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and X2 translated to the notation in this paper, the RW metric has the time 
measure dti and the length measure dx\ = a{t)dx and the transformed time 
measure dt2 = dT and a length measure dx2 = dR. Then the criteria for 
clock physicality becomes dit2 = Tt = (?2ti = = l"^ /Tt (i.e., ^ = 1) and 
for ruler physicality 81X2 = Rx/a{t) = 82X1 = a{t)xR = J^a{t)/Ry^ (i.e., 
B = 1). Thus, the criteria for physicality at any distance are 

rt = 7,i?^ = 7a(t). (26) 

If we expand Tt and in a power series in %, we will require that the 
criteria in eq 26 are valid to the lowest non-zero order of %• This makes 
quantitative the meaning of "close to" as compared to "at" the origin. We 
will find later that there are an infinite number of coordinates systems which 
meet this requirement. 

Since V will vanish at the origin, we will require the coordinates to be 
physical to the lowest order of V and/or 

Tt = l + \v\ (27) 

and 

R^ = a{l + \v^). (28) 
if we integrate eq 28 we obtain 

1 

R = ax + -j^ aV^8xt, (29) 

where the symbol dxt signifies that the integral is to be carried out with 
respect to x with t held constant. We can neglect the integral in eq 29 since 
it is of the order of x^. The resultant distance is the proper distance to 
which all measurements of distance reduce close to the origin [4]. 
We can find V from eqs 14 and 29: 

V = -dx + 0{x^), (30) 

where the dot is the derivative with respect to t. Thus V vanishes as the 
first power of x with no second order term. This result can also be obtained 
without using eq 22 by invoking the Bernal criteria directly on eqs 9 and 16, 
then expanding D in powers of V and cT^/aTt, since both must vanish at 
the origin when C = 0, and then finding V from eqs 14 and 29. Similarly, 
we can find an expression for T, using eq 22: 

T = t+ T^dxt = t- -VTtdxt. (31) 
Jo Jo c 
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If we partially differentiate by t, we get another expression for T^, using eq 
30: 

If we equate the two expressions for Tt (eqs 27 and 32), and neglect a higher 
order term, we obtain an integral equation for c: 

To get an explicit expression for c, partially differentiate both sides by %: 

«=X = X^(-). (34) 

Cancel multiply by a, change the variable dt to da/ca, and integrate to 
obtain 

c/co = = a-E, (35) 
do 

where a is the normalized scale factor 

a = a/ao, (36) 
and E is the normalized Hubble ratio H{t) 

£iQ tiQ a 

The subscript denotes the value at t = toj the present time. We will take 
Co to be unity, so that c{t) will be measured in units of cq. The field equation 
(sect 5) will enable us to evaluate a and E and thus c{t). 

Thus, we have shown that if transformed physical coordinates exist, they 
must have the variable c{t) shown in eq 35. In Sect. 4 we show that there 
exist an infinite number of such transforms. 

The physical coordinates have a generalized Minkowski metric close to 
the spatial origin (C = 0, A = 1, B = 1, ar ^ ax — > R)'- 

ds^ = dT^ -dB? - B?dio'^, (38) 

This metric is invariant under a generalized Lorentz transformation between 
a frame [dT' , dR') moving at a radial velocity v relative to a second frame 
{dT,dR): 

df' = ^'it)idf-^dR), 

}' (39) 
dR' = i{t)i-^dT + dR), 
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where j'{t) = 1/^/1 - /c{ty. Of course, because of spherical symmetry, 
this applies to a frame moving in any direction for differential distances 
measured in the direction of the velocity v. The generalized Minkowki metric 
and generalized Lorentz transform become the usual ones for events that 
occur in times short compared to the cosmic time t. This, of course, is a 
requirement for physicality, since it so basic to all of our physical laws. We 
could have used the generalized Minkowski metric as an alternate physicality 
requirement to get the same c{t). 

If we had used dr = ds/ c{t) as the invariant line element to define the 
covariant transform, close to the origin we would have obtained 

dT^ = dT^ - -^{dR^ - R^du'^), (40) 
c{ty 

with the same c{t) as above. But the time dependent coeficients of the space 
elements would provide a major complication if we tried to use dr in the 
field equation, which cannot be simply fixed by incorporating it in dR and 
R. This is a good reason not to use it. 

It should be observed that this variable photon velocity was derived with- 
out the use of the field equation. The existence of the variability is caused 
by the homogeneous and isotropic symmetry assumed for the universe, the 
validity of the covariant transform (including the invariance of the line ele- 
ment ds) , and by the concept of physicality extending to the lowest non-zero 
order of x, but not by the gravitational energy density, whose only effect on 
c{t) is how it affects the rate of change of the cosmic scale factor a. 

4 Complete physical transformed coordinates ex- 
ist 

We show in Appendices A and B that there exist an infinite number of 
coordinate systems which satisfy our physicality requirements. Appendix 
A derives the equations for diagonal transforms using phyical clocks for all 
distances {A = 1) and Appendix B for physical rulers {B = 1). Because 
the two sets are different from each other, they show, as we would expect, 
that it is not possible to have diagonal physical rulers and physical clocks 
simultaneously for all values of t,x (except for an empty universe). 

I have been able to find analytic solutions for all the physical clock trans- 
forms, but not for all the physical ruler transforms. However, because there 
is an interest in having a physical description for distance in the universe, we 
have a special interest in the physical ruler transforms. I have been able to 
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find a similarity solution for physical rulers for the special case of Euclidean 

space with $7 = 1. The result is displayed in Fig. 1. Here we have used the 
field equations with the generalized time (Sect 5) to derive the equations 
for a{t) and c(t). Plotted against the time at the origin (cosmic time t) are 
galaxy distance (constant x) for various galaxies and the photon distance 
for light emitted by the galaxy % at t = and reaching the origin at t = ti. 
Notice that there is a horizon at x = 4 for the visible universe ti = to, and 
that light comes monotonically towards us from all galactic points. This 
photon path has a slope of c = (to/^i)''^^^ close to the origin where the dis- 
tance R and time T = t are both physical. Although the distance uses 
physical rulers, the coordinate system as a whole may not be physical for 
> 1/4, i? > t^/^, shown by the dashed line in Fig. 1. This excludes 
much of the paths at early times. 



5 Revised Equivalency Principle for a variable pho- 
ton velocity. 

To accomodate c(t), it would seem appropriate to revise the Equivalency 
Principle as stated by Misner et al[6] to read "to any and every local gen- 
eralized Lorentz frame, anywhere and any time in the universe, all the non- 
gravitational laws of physics must take on their generalized special relativis- 
tic forms", viz. generalized Lorentz transforms (eq 39) with the generalized 
Minkowski metric (eq 38). The use of t or T in the field equation accom- 
plishes this: 

StvG 

G^,^ + A9^^ = -^T^u, (41) 

where G^i, is the Einstein contracted curvature tensor determined from the 
coefficients g^i, multiplying the differentials in the metric of eq 4, or of eq 8, 
A is the cosmological "constant" possibly representing some kind of vacuum 
energy density, G is the gravitational "constant" and T^i, is the energy- 
momentum tensor. Locally, this will give the generalized Minkowski metric 
(eq 38) and generalized Lorentz transform (eq 39) for dT and dR close to 
the origin, as well as for dt and adx- 

It will be convenient to use the variables t, r, 6, cf) in the generalized RW 
metric (eq 4): 

^ ^£2 _ a\^^ + rHe^ + r^cos'edcp'). (42) 
1 — Kr^ 

so the g^u of eq 41 are the coefficients of eq 42, exactly the same as the usual 
RW metric with constant c = 1 in the variables t, r, 9, (p. We will make the 
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usual assumption that the universe is an ideal fluid with an energy density 
of p(? and pressure p, so that we can write the two significant field equations 
[6, page 729] for a{i) as 

3d2 3A: , SttG 

— + ^ - A = -^p, 43 

and 

'd k , SttG , , 

+2- + ^ + ^ - A = ^p, 44 

where the dots represent derivatives with respect to t. 

Now, the first equation can be multiplied by a^/3, differentiated, and 
subtracted from da^ times the second to give 



d ,Gp 3G , . 

^(io'l—^-^f-i^^A, (45) 



where we have allowed for the possibility that G and A may be functions of 
c{t) . For small p and A, 

Gpa^/c^ = constant (46) 

If the energy density consists of n particles per unit volume of mass m, so 
p = nm, then the conservation of particles requires na^ be constant. This 
makes 

Gm . , 

— ^ = constant. (47) 

Even when p and A are not negligible as in the early universe, it is reason- 
able to assume that matter particles are conserved well into the relativistic 
particle era and that the relation of constants in eq 47 will continue, so that 
for matter density eq 47 remains valid. 

6 Calculation of a{t) and c{t). 
Following Peebles [7, page 312], we define 

3Hi4 

and 



^ - ^ (48) 



Qr — tt2„2 (49) 
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and 

SlA ^ ^. (50) 

For very small a there will also be a radiation term whose value will have 
to be reexamined for its dependence on c{t). For this paper it will not be 
included. 

The Qs are defined so that 

n + nr + ^A = 'i- (51) 
Then the normalized Hubble ratio E in eq 37 is determined by eq 43: 



which allows us to evaluate c{t) = aE. Ai t = to: oq = 1; E{ao) = 1, and 
Co = 1. 

The physical time t becomes 



coHot 



f" da da. 

k ^^Jo ^^^^ 



where co = 1 is written explicitly for clarity. For = 1: c = a ^^'^ = 
{to/ty^^, and coHot = whence coi^oto = 1/2, a = (t/to)^/^- 

For experiments attempting to measure the variation of the photon ve- 
locity at the origin at the present time, the derivative of c{t) (eq 35) will be 
more useful: 

-^\-^] =i-ln-nr = ~ + n^. (54) 

coi?o LcaiJt=to 2 2 

Notice that this fraction is negative when matter dominates, and goes from 
zero at zero density to —1/2 at the critical universe density. A vacuum 
energy density opposes the gravitational efi'ect of matter, and as the vacuum 
energy density increases, the slope goes to zero before reversing from a 
decreasing function of time to an increasing function. 

7 Other physical constants. 

Many physical quantities depend on c. The fine structure constant a j in SI 
units [9] is 



e2 



(55) 
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and the Rydberg constant expressed as a frequency is 



RooC = oif—^ = „ ■ (56) 

Note that the 47reo is often omitted in the fine structure constant since it is 
unity in Gaussian coordinates, but it is essential here if we are to consider 
a variable c{t). 

Atomic clocks have reached a stability and precision to demonstrate that 
a/ and i?oo remain essentially constant with time [10] [11]. Prom this we 
can conclude that 

^ oc c(t), (57) 

and 

Another relation we can count on comes from electromagnetic theory: 

eo/xo = (59) 

Lastly, if the use of generalized cosmic time in the field theory is correct, 
we obtain (eq47): 

GmQcc{tf. (60) 

If e and h are constant, we could conclude that = l/c{t), hq = l/c(t), 
m oc l/c(t)^, and G oc c(i)^, but this is only one possibility. The actual 
variation will depend on comparison with experimental data, such as super- 
novae la, cosmic background radiation, gravitational lensing, and dynamical 
estimates of galactic cluster masses. For this, of course, we need a theory of 
the c{t) dependence of each of these. 



8 Conclusions 

I have investigated what coordinate systems, other than the Robertson- 
Walker (RW) system, there might be to describe a homogeneous and isotropic 
universe that are physical close to the origin, i.e., are derived by a covariant 
transform from the RW coordinates and use physical clocks and rulers. Ho- 
mogeneity requires these coordinates to have diagonal metrics. I have shown 
that all such systems must have the same variable photon velocity close to 
their origin, which is consistant with the RW system also having the same 
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variable photon velocity. Since this is determined by physical coordinates, 
it is reasonable to suppose that this represents the physical photon velocity. 
I find that the analytic expression for it is c{t) = aH, the cosmic scale factor 
times the Hubble ratio in appropriate units. For a homogeneous universe, 
the origin can be placed at any point in the universe, so this variable photon 
velocity should be applicable to all physical processes everywhere. 

The existence of the variability is caused by the homogeneous and isotropic 
symmetry assumed for the universe, the validity of the covariant transform, 
and by the concept of physicality of the transformed coordinates close to 
the origin. If we deny the variability of the photon velocity, we must find a 
flaw in one of these concepts. It is not caused by the gravitational energy, 
whose only effect on c{t) is how it affects the rate of change of the expansion 
coefficient a. 

A variable c(t) no longer keeps the local Lorentz transform invariant. 
However, we can make it invariant by using a generalized time for which 
the differential is c(t)dt. We can then create a covariant derivative and 
a local generalized Minkowski metric in order to use Riemann algebra in 
Einstein's field eqation. With this I find that the Hubble ratio defined with 
this generalized time has the same dependence on the cosmic scale factor as 
it does for a constant photon velocity. 

Some other physical "constants" must also be changing with cosmic time, 
but its hard to say which ones. I have suggested some constraints on this 
variability, including keeping constant the fine structure and Rydberg con- 
stants. 



Appendix A: Transformed coordinates with physi- 
cal clocks 

For diagonal coordinates with physical clocks at all t and x, ^ = 1. Thus, 
eq 23 becomes 

Tt = 7- (61) 

We proceed by finding a differential equation with only V as the de- 
pendent variable. Thus, we can write a formula for T, using this and eq 
22: 

T = t+ r T^dxt = t + r{--lV)dxu (62) 
Jo Jo c 

where we have used the boundary condition that at % = 0, T = It can be 
partially differentiated with respect to t (giving 7) and then with respect to 
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X and with the use of eq 10, noting that d'y = ^^VdV and 1 + F^7^ = 7^, 
we obtain a PDE for V: 



R 



R 

This can be rewritten as 

dVR _ d{%)R 

F(i-y2) f 



(64) 



where the subscript on the partial differential indicates the variable to be 
held constant. This can be integrated with an integration constant In/t. 
Since the integration is done at constant i?, then k = k(-R), and inversely, 
R = R{k). Integrating eq 64, we get 

V = - , (65) 

where the sign of k will be positive for an expanding universe, where the x 
points will stream out radialy past a point at R. 

At this point, R is an unknown function of k. The various possible 

coordinate systems which solve our PDEs are characterized, in large part, 
by the function R{k). But for all, in order for V to vanish when R = 0, k 
must also; so always 

k(0) = 0. (66) 

One observation we can make without any more work is that, as long as 

k{R) remains finite, V goes to —1, and V goes to — c(t), for d{t) = 0, i.e. 
for t = 0, the horizon. This is different from any of the three definitions 
of distance assuming a constant photon velocity at the origin treated by 
Fletcher (including the proper distance) [5], for all of which V became — 1 
at a finite t instead of at t = 0. 

Let us now look at lines of constant k,{R), i.e. constant R, in t,x space. 
Eq 10 can be integrated for x with use of eq 65 at constant k to give the 
following: 

x{t,K)=r ^f^^- . (67) 



For an open universe, we have set the upper limit at 00, becaause we expect 
that if R is kept constant the galactic point x that will be passing any given 
R will eventually approach zero as RW time t approaches infinity. This 
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limit may be valid even for a closed universe, since c{t) goes to zero for long 
times, making the time go to infinity perhaps before the first maximum of 
a periodic R. 

At this point, we have obtained V = V*(t, k) from eq 65 and have also 
obtained the function k). We can in principle invert eq 67 to obtain 
K in terms of t and x- = K(t,x)- This gives us the velocity function 
V{t,x) = y*{t, K{t,x)- If the function R{k) were known, we would then 
also have R{t,x) = R{K(t,x))- 

T{t,x) can be found by noting from eqs 22 that 



aV 



-7 



K. 



(68) 



By substituting eq 68 into eq 62, and integrating over k instead of x by 
dividing the integrand by the partial of eq 67 with respect to k, we find an 
expression for T{t,x)' 



T{t,x)=t + 



1 



^1 + "^""^ 



(69) 



where k is put equal to K{t, x) after integration at constant n in order to 

get T{t,x). 

This completes the solution. Since k{R) can be any funtion that van- 
ishes at the origin, there thus exist an infinite number of solutions for our 
transformed coordinates with A = 1,C = 0. 

To determine physicality, we will next find a/R^ close to the origin. 
R^ = k^/k'{R) can be written in an inverted form by taking the derivative 



of eq 67 with respect to k: 

-1 



1 



= K 



f dK 



dx 



t a2(l + ^)3/2 



(70) 



To be physical ^a/Ry. = 1 as R approaches 0. Putting 7 = 1 and k{0) = 
k'{0)R{0) = in eq 70, multiplying by a, setting the resultant equal to one, 
and changing the integration variable from t to a{t) gives 



i=„'(o)rf|i. 



(71) 



remembering that the dot indicates differentiation by t. k'(0) is a constant 
to be determined by c{to) = 1. Note that the integral of eq 71 is independent 
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of the functional form of k{R), and is therefor the same for all k{R). In fact 
it is an integral equation for c{t) . By differentiation of both sides of eq 71 
by a, we can obtain 

c(t) = (72) 

which, as we should expect, is the same c{t) (see eq 35) we showed for all 
physical coordinate systems for k'(0) = clq. 

It is interesting to consider the limiting case of a zero density universe: 
$7 = 0, = 1, o.qHq = 1 (eq 49). Eq 51 makes a = Hq for all t,x- 
Integrating gives a = t. Eq 35 makes c = 1. Eq 67 gives x = csch~^t/K, or 
K = K(t, x) = t sinhx- We can then find from eq 65 that V(t, x) = — tanhx 
and from eq 11 that 7 = coshx so that 

dRKy (IR 

c= — = ^ ^ = ^ = 1- (73) 

7a an 'ft an 

Thus the physicality condition is met for all R with R = K and A = \,B = \, 
so that the complete transform with eq 69 becomes 

R = tsmh.x, 

T = tcoshx- ^ ' 

These coordinates have been known ever since Robertson [12] showed that 
this transformation from the RW co-moving coordinates at zero density 
obeyed the Minkowski metric. What is new is that this solution was derived 
from the equations we obtained for our physical clock transforms with A = \. 
It can also be obtained from the physical ruler transforms {B = 1) since eqs 
63 and 79 for V become identical with Vt = and a/c = a = t. 

Appendix B: Transformed coordinates with physi- 
cal rulers 

For diagonal coordinates with physical rulers for all t and x, S = 1, that is 

Rx = «7 (75) 

By integration we find 

f'X 

R = a jdxt, (76) 







and partial differentiation with respect to t gives 



rx rx 
Rt = a -fdxt + a ltdxt- (77) 

JQ Jo 
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We can then find V from eq 14 as 

cry C7 



rx rx 

a / -ldxt + a I ltdxt 
. Jo Jo 



(78) 



This is an integral equation for V. It can be converted into a partial differ- 
ential equation by multiplying both sides by 7 and partial differentiating by 



7 



+ -VVt 
c 



1 da 
c dt 



(79) 



Note that this is substantially different from the eq 63 for V that we obtained 
for physical clocks. This means that it is not possible to find diagonal 
transforms with both physical clocks and physical rulers for all values of 
t and X (except for 17 = 0). It is possible to have either one or other be 
physical at all t and x with the other being physical only close to the origin. 
I have not found it possible to obtain explicit general integrations of eq 79 as 
we did for physical clocks. However, I have found a similarity integration for 
the special case of = 1 where a = ao{t/tQy^'^ and c = (to/t)^^^. To simplify 
notation let us normalize time to t/to — > t, a/ao —>■ a, and xco^o/oo x, 
T/to T, R/coto R, and let W = -V. Eq 79 then becomes 



(80) 



This can be converted into an ordinary differential equation (ODE) by let- 
ting 

n = x/tV4 (81) 

so that eq 80 becomes 



W'{l + '^) = \{l-W% 

where the prime denotes differentiation by u. 

Similarly we can find ODE's for T and R by defining: 



(82) 



T = tq{u), 



(83) 



and 



R = t^/^s{u), (84) 
where q{u) and s{u), from eqs 22 and 14, are given by the coupled ODE's: 

uW. 



qW, 



(85) 
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and 

s'{W + -) = -s (86) 

For small values oi u, W = u/2, g = 1 + s = u, and i? = t^^'^x = ax. 

We also note that Tt = 1 + confirming that these coordinates have a 

physical clock close to the origin, justfying c(t) = t~^^^. 

An alternate approach would be to start with c{t) unknown, but of the 
form c = t~^. This will also yield a similarity integration of eq 79 with 
Tt = 1 + 2bW^ for small u. In order to make the clocks physical close to the 
origin, Tf has to equal 1 + so that 6 = 1/4, confirming that for these 

diagonal coordinates to be physical at the origin (x da/ dt, oi c oc da/dt. 

As t ^ 0, u ^ oo, and W'^ ^ 1 - 700/{u + 4)^, q .0014(m + 4)^ and 
s .012('u+4)3. T and R both remain finite at this limit with T .0014a;^, 
and R .012x^. The fact that T does not go to zero when t goes to zero 
results from synchronizing T with t at the origin, so that at finite R, where 
clocks run slower than t, they have different starting points. The distance R 
vs the time t of the observer at the origin can be found from the numerical 
integration of the coupled ODE's, eqs 82, 84, and 86. The paths of galactic 
points are those for constant x. The path photons have taken reaching the 
origin at ti is found by calculating Xp vs t and using the tranform to R. 
Thus, for = 1 

x„ = /"* -dt = 4(t:;/^ - tV4) (87) 

Galactic and photon paths arc shown in Fig. 1. Although the distance uses 
physical rulers, the coordinate system as a whole may not be physical for 
W^/2 > 1/4, n > .90, R > Mt^/^. This is shown by the dashed line in Fig. 
1. This excludes as non-physical much of the paths at early times. 

There are an infinite number of similarity solution for the physical rulers, 
obtained by different ratios of x to x, which are still physical close to the 
origin. These give different numerical limits of T and R at small t, and are 
thus different functions of t and x- There are also a corresponding number 
of similarity solutions for physical clocks. I have numerically integrated one 
of each with twice the ratio of x to a; of that derived above. All three are 
very similar in both galaxy paths and photon paths for the same values of 

X- 
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Figure 1: Physical ruler distances for r2 = 1 
Distance R in units of coto for various galaxy paths (constant x = cotox/o-o) 
and for various light paths plotted against the normalized time on clocks on 
earth {t/to). The visible horizon is the path that the farthest galactic matter 
{x = 4) takes that is visible today (t = Iq). The light paths labeled with 
X are the paths which the photons take after emission by the galaxy x at 
i = 0. Notice that the slope of the light path close to t = to is cq = 1, where 
the coordinates R and T = t are both physical. For light paths arriving 
at some earlier time ti, the slope close to i? = is c(ti) = (ioAi)"*"^^- 
The dashed line in Fig. 1 shows the upper limit of assumed physicality 
y2/2 < 1/4, i? < .9^3/4. 
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